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Abstract 

We study the (tree- level) weak radiative decays of B mesons. We present a numerical 
estimate for the inclusive b — ► ^7(7) modes based on the free-quark decay. We then 
review what is known for the B — > D*j modes in the framework of heavy quark effec- 
tive theory and chiral perturbation theory. Finally, we extend these ideas to the double 
radiative decay modes B — > D77. We find that the b — ► X c 77 rate is about an order 
of magnitude larger than the corresponding b — ► X s j-f rate. We also find the branching 
ratio for the B — > D77 mode with most favorable CKM angles at the few xlO -8 level, 
comparable to predictions for B — > K'yy. 

PACS numbers: 12.39.Hg, 13.25.Hw, 13.40.Hg 



*Elcctronic address: oaanti02@iastatc.edu; valencia@iastate.edu 



1 



I. INTRODUCTION 



The radiative penguin b decay of the form b — > 57 has received much attention in 
the literature because it is sensitive to certain types of physics beyond the standard 
model. The HFAG quotes an average for the measured branching ratio B(b — > 
57) = (354+2s) x 10~ 6 (H- Recently the double radiative decay mode b — > 577 
has received some attention in connection with the possibility of measuring it at a 
Super-B factory at the 10 -7 level 0]. This is the level at which it is expected to 
occur in the standard model 0, 0|- 

There is another type of radiative b decay modes in which the weak decay pro- 
ceeds via the charged current at tree-level, b —> X^^) and b — > ^7(7) . These 
modes are expected to be dominated by standard model physics and have received 
much less attention. We begin this paper with a simple numerical estimate for these 
modes suggesting that b — > X c yy can also be observed at a Super-B factory. We 
then turn our attention to some of the exclusive modes. 

The charged mode with one photon, B + — > -D + *7, has been studied before in 
the context of heavy quark effective theory as a potential probe for V u b by Grinstein 
and Lebed 0. The neutral mode with one photon, B° — > -D°*7, has been studied 
by Cheng et. al. in a slightly different context. We first review what is known 
about these modes using the framework of heavy quark effective theory and chiral 
perturbation theory and present new results that include the effect of intermediate 
positive parity states. We find that these intermediate states have an important 
effect on the overall rates, particularly for the charged modes, due to a partial 
cancellation that occurs in the leading order amplitude. 

We then extend these results to the case of the double radiative decay modes 
B — > D77. We include the lowest lying positive parity states in our calculation and 
find that they play an important role, just as they do in the single radiative decay 
modes. The double radiative decay modes present, in principle, the opportunity 
to study the heavy quark and chiral expansions by looking at different kinematic 
regions. In practice, however, this may not be possible due to the small rates. We 
find that the mode with the most favorable CKM angles can occur at the few xlO -8 
level and is perhaps observable at a Super-B factory. 



II. b -» X c7 AND b -> X c77 

We begin by discussing the inclusive radiative decays. We work in the free 
quark approximation where b — > X c y and b — > X c r y r j arise from the tree level quark 
processes b — > cudy and b —> cud^yy (we take V u a = 1). Our goal is to compare these 
processes to the one-loop processes b — > X s y and b — > X s yy that have been studied 
at length in the literature. We will limit ourselves to a numerical estimate with the 
aid of CompHEP 0. 

The process b — > X^ is to be compared with the inclusive b — > Xg^. Accordingly 
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we impose a cut on the photon energy 1.8 GeV < E 1 < 2.8 GeV corresponding to 
the energy range studied by the Belle collaboration We also impose a separation 
cut that requires a minimum angle Q m i n between the photon and the final state 
quarks. With values of 5° < Q m %n < 20° for this minimum angle we obtain rates 
r(b — > cud'j) between 3 x 10~ 18 GeV and 5 x 10 -18 GeV. To estimate the branching 
ratio for the inclusive process we then use 

B(b -> X c >y) = — ^B(B + -> X c e + is) exp 

T(b -> ce u) th 

~ (7-11) x 10~ 6 . (1) 

To obtain this number we used V c b = 0.0413 jsj, B(B + — > X c e + u) exp = (11.15 ± 
0.26 ± 0.41)% [ltj and quark mass values m b = 4.8 GeV, m c = 1.5 GeV which 
are the ones used in the theory estimates of b — > X s / -f-f. We show in Figure ^ the 
photon energy spectrum for the case oid m i n = 5°. The characteristic bremsstrahlung 
spectrum falls rapidly with the photon energy and for this reason B(b — >• X c ^) is 
much smaller than the penguin mode B(b — > X s ^) = (3.3 ± 0.4) x 10~ 4 in this 
energy range. 
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FIG. 1: Photon energy spectrum in b — > X<pf with 6 m i n = 5°. 

Next we consider the double radiative mode B — > X c / -f~f in a similar way. We 
wish to compare it to the penguin process B — > X s 77 which has not been observed. 
We use instead the theoretical calculation of Reina et. al. 0,0] which finds B(B — > 
^77) ~ 3.7 x 10~ 7 (including LO QCD corrections and using the quark masses 
mentioned above). They also require that the energy of each photon be larger than 
0.1 GeV, that the photon pair invariant mass be larger than 0.1m&, and that the 
photons be separated from each other and from final state quarks by at least 20°. 
They find a spectrum that is sharply peaked at low M 77 invariant mass. 
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For our calculation we estimate B — > X c 77 from the leading tree-level contribu- 
tion b — > cud'yy. With the same cuts used by Ref. 0, 0] for B — ► X s 77 we obtain 
r(6 -> CMc/77) = 1.9 x 10~ 18 GeV and this goes up to 3.2 x 10~ 18 GeV if the angular 
cuts are relaxed to 10°. Proceeding as in Eq. ([IJ we arrive at 



B(B -> X c77 ) ~ (4.2 - 7.2) x 1(T 6 . 
In Figure 121 we show the two photon invariant mass distribution. 



(2) 




FIG. 2: Mr distribution for the process 6 — ► X c 77 with the cuts described in the text. 



Once again we see a spectrum that is strongly peaked at low invariant mass. 
Comparing our results to those of Ref. 0, 0] it is clear that B — > X c 77 can be a 
significant background to B — > Xg'-f-f. Notice that the total rate for the double 
radiative decay mode is comparable to that for b —>■ X c ^. This of course is due to 
the much tighter photon energy cuts imposed in that case. 



III. EXCLUSIVE MODES AND HQET 



In this section we collect the ingredients necessary to calculate the amplitudes 
for the exclusive modes guided by heavy quark effective theory, combined with SU(3) 
chiral perturbation theory . 

1 Similar methods were used in Ref.^lJ to calculate rare D° decay to two photons 



4 



A. Radiative decays of heavy mesons 



The strong interactions involving the heavy meson (0~,1~) doublet (the light- 
est pseudoscalar and vector mesons containing one heavy quark) and light pseudo- 
scalars are described by the effective Lagrangian (we drop the subscript v that 
indicates the velocity of the heavy meson to simplify the notation) 

£ = -iTr(HPv ■ D ba Hf ] ) + gTriH^Hf^^A). (3) 

In Eq. © we use the standard notation in which 

• The heavy pseudoscalar and vector meson fields with heavy quark Q and light 
anti-quark q a are destroyed and created by the field HP and its Hermitian 
conjuagte. They are given by 

H (Q) = l+t (P « V _ P (Q) T5 ) (4 ) 

which transforms under chiral symmetry as — ► H^U^. The Hermitian 
conjugate matrix = 7 i?a Q ^7o- 

• Mesons containing heavy anti-quarks Q and light quarks q a are destroyed and 
created by the fields 



HP = (Pi? V " Pi 0) 75) 



1 - f 



SP = loHpho (5) 



• The pseudo-Goldstone boson octet, <p is incorporated into a 3 x 3 unitary 
matrix £ = exp (2i 4>/f n ) which transforms under chiral symmetry as £ — > 
LluR) . In Eq. (JHJ) they enter through the matrix £ where £ = £ 2 and the 
transformation properties of £ under chiral symmetry £ — ► L£W = U£R> 
define the matrix U. The charges of the light quarks appear through the 
diagonal matrix Q with entries 2/3, —1/3, —1/3. 

• The chiral covariant derivative and axial current in Eq. J3J) are given by 

D» ab = 8 ab d» - V£ = 8 ab d» - + £d^) ab 

< b = ^d^-td^) ab (6) 

The leading order electromagnetic coupling is obtained from Eq. (JHJ) by minimal 
substitution. However, the couplings arising from this procedure (for charged B 
and D mesons) do not contribute to the processes B — > or B — > .D77 as can be 
seen by explicit computation. The lowest order electromagnetic coupling that will 
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contribute to these processes is the transition magnetic moment 0,0000. 
For mesons containing a heavy quark it can be written as, 

C em = -^^TriH^a^HP)^ ~ ^Tr{H^ H { b Q) a, u <£ ba F^), (7) 

where = \{gQ£ + £Q&). 

The coupling consists of two terms corresponding to the decomposition of the 
electromagnetic current into heavy and light quark parts, 

Ha = e Q ^ + eafi® = ^ + $ = + (8) 



Ao A 



a 



where e a is the charge of the light-quark. The heav y q uark contribution at leading 
order in the 1 / mq expansion is given by Aq = mq |l9l . 0] • The lig ht quark contri- 
bution in the SU(3) symmetry limit is usually called Aj* = j3 |l4j . This constant 



has been estimated in vector meson dominance models [21| as well as in chiral quark 
models 0] (along with SU(3) breaking corrections). The leading SU(3) violations 
have also been calculated |l6[. More generally, the heavy quark contribution when 
the mesons have different velocity is modified to fii^ = eg((w)/mQ 0,113 where 
uj = v ■ v' and is the Isgur and Wise function. 

Eq. (0) generates the following amplitudes, 

M(D*( Vl ) - D*(7fe)7(g, e)) = (q ■ Vl e* ■ rf 2 - q ■ r&e* ■ V i) , (9) 

where we have defined /j,d* = (l^a — Ha^)- Analogous relations with /j>d,d* 
Hb,b* then hold for the B system. In the heavy quark and SU(3) limits, the magnetic 
moments consist only of the light quark contribution given by \i a = e a (3. For our 
numerical estimates we will use the leading magnetic moments as well as the mag- 
netic moments in three models tabulated in Ref. 0]: "xLM" a chiral loop model; 
"VMD" a vector meson dominance model; and "RQM" a relativistic quark model. 

The magnetic moments from Eq. (J7J) are defined for on-shell transitions between 
a vector and a pseudo-scalar of the same mass. In our calculations one of the mesons 
will be off its mass shell and the corresponding form-factors should be evaluated at 
k 2 ~ —5m 2 for the single radiative decay modes: e/i = «7a/(0) — > e/i(?A/(— 5m 2 ) where 
5m = mb — m c . On general grounds one expects the form factor to change over a 
characteristic scale Aqcd, so there is large uncertainty associated with the use of 
the on-shell form factors. Formally the results we obtain, correspond to the so called 
generally low velocity (GL) limit in which j^| 

5m ~ Aqcd < m b . (10) 

Alternatively one can model the momentum dependence of the form factors as was 
done in Ref. for the neutral modes. We will not include a momentum dependence 
of the form factors in our estimates but instead introduce two additional effects. 
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First we will keep certain terms that are formally of order 5m /nib arising from spin 
one propagators as described later on. We will also consider additional heavy meson 
intermediate states; the positive parity P-waves of the system Qq. However, we will 
neglect higher total spin resonances. 

The positive parity states that we include are predicted by HQET to lie in 
two distinct multiplets: (0 + , l + ) = (i?o, B\) and (1 + , 2 + )=(Bi, B2). Generically we 
will refer to them as B** and will discuss the case of B mesons for definiteness, 
with corresponding results for D mesons also used in our calculation. The velocity 



dependent fields are introduced in a manner similar to the field H in Eq. (J3J) 24 , 



S = i(i + fl# l75 _fl 
T» = -(1 + 
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The electromagnetic transitions between a member of these doublets and a member 
of the (CT, l - ) doublet have also been discussed in the literature. Analogously to 
Eq. (J7J) the couplings receive contributions from the heavy and light quark currents. 

The heavy quark contribution to the (1 + , 2 + ) to (0~, 1~) transition in the charm 
case can be found in [2^. A simple way to reproduce those results consists of using 
the matrix elements for < B**\V il \B > obtained by Isgur and Wise |23| to match 
the effective Lagrangian £ ~ Tr [H^T^ QijF a p). Similarly one can start from the 
Isgur and Wise results and impose gauge invariance for the kinematic conditions we 
consider to obtain the heavy quark contribution to the (0 + , 1 + ) to (CT, l - ) transition. 

The light quark contributions can be determined at leading order in chiral per- 
turbation theory in terms of two unknown constants, 

C = ~^Tr {HiT?^ QijF a p) - ~^Tr (H^a^Q^) + h. c. (12) 

A 3/2 4A l/2 

Combining the heavy and light quark contributions for the T doublet leads to 
the amplitudes (we include only those that do not vanish at leading order in the 
1/tuq and chiral expansions) 

M(B 2 (v,e 2 )^B*(v,e v )y(q,e)) = -2V3efi T B (ef q a e* vp e* ■ v - ef e^e*g • «) 

^( J B 1 (t;,e 1 )^ J B*(t;,e y ) 7 (g,e)) = i^le^q^e^ 

M(Bi(v, ei ) -> B(v)j(q, e)) = 2y/2 ef i% (q • e lV ■ e* - q ■ ve l ■ e*) . (13) 

Similarly for the S doublet we obtain 

M(B 1 (v,e 1 ) ^ B^v^vMq^)) = -ie^e^q^e^p 

M(Bi(v, ex) -> B(v)j(q, e)) = e/x| (q ■ e x v ■ e* - q ■ vt x • e*) 
M(B (v) ^ B*(v,e v )y(q,e)) = e^ s B (q ■ e\v ■ e* - q ■ ve* v ■ e*) . (14) 
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Corresponding expressions for charm are obtained with the obvious replacements 
/i^' 5 — > ^ S ■ For bottom these effective coupling constants are 

-H^+i)- (i5) 

For the light quark contributions, A' 3 ^ 2 corresponds to y/3/ (/ — /') of Ref. j2^|. Using 
their estimate for the D system, we take A' 3 , 2 ~ (2.75 — 3.5) GeV. Similarly, A^ 2 
corresponds to A' of Ref. |27j where it is estimated that A' x , 2 ~ 1.25 GeV. For the 
heavy quar k contributions we have used the Isgur-Wise functions r 1 / 2,3 / 2 estimated 
in Ref. [3 to be r 1 / 2 ' 3 / 2 (l) ~ 0.24 (£ 3 / 2 of Ref. Q corresponds to \/3r 3 / 2 ) . More 
recent estimates of these functions using the light-front formalism [22| and lattice 
calculations indicate larger values, particularly for r 3 / 2 (l) up to 0.61. These 
estimates also indicate that these functions vary by less than factors of two over the 
range of u that is kinematically allowed in B — > D77. 

Finally we will also need (0 + , 1+) to (0 + , 1+) and (0 + , 1+) to (1+, 2+) electromag- 
netic transitions. The former receives heavy and light quark contributions, whereas 
the latter only receives light quark contributions at order 1/mq. They can be ob- 
tained from the Lagrangian 

C = -^Tr{S a a, u S a )F^-^Tr{S a S b a, v Q ab Fn-^Tr (S^Q^) . 
4m Q 4A 1/2 M A 3/2 3 

(16) 

The new constants Ai/2,3/2 are not known. For our numerical estimates we use 
Ax/2 ~ A a and A3/2 ~ A3/2. With this choice the magnetic transitions within the S 
multiplet are the same as those between members of the H multiplet, for charm for 
example one has 

M(Di(rj) —> D j(q,e)) = -ie^e^^q^ 
M(Di(r]i) -> D 1 (r) 2 )j(q, e)) = (q ■ r^e* ■ rf 2 - q • r/*e* • 771) , (17) 

and corresponding expressions for bottom. The (0 + , 1 + ) to (l + ,2 + ) vertices for 
bottom are 

M(B 2 (v, e 2 ) - B 1 (v, e v )>y{q, e)) = i2ef B s [f^efqpv^ + ef q^e*^ 
M(Bi(v, e<) -> Bx(y, e / ) 7 ( ? , e)) = -e/i% s ^ (q ■ e,e* ■ e} - q ■ e)e* ■ e*) 



M(B 1 (v,e 1 )^B (v) 1 (q,e)) = i2ef B s ^ -e^q^^, 
and corresponding expressions for charm. We have defined 



H™ = (19) 

A 3/2 



s 



B. Weak transitions 



Within the standard model, the effective weak Hamiltonian responsible for the 
AB = 1 transitions at tree-level is (with di = d or s) 

n = ^= [v; b v udi {cm + c 2 Q n 2 ) + v: b v cdi (c.qi + c 2 q c 2 )\ + h. c . (20) 

where the neutral modes B° — > D°(D *) are mediated by Q™ = (bdi)v-A(uc)v-A 
and Q 2 = (bc)v-A(udi)v-A', and the charged modes B + — > D + (D + *) are mediated 
by Qt = (J>d>i)v ~a(cu)v -a and Q 2 = (bu)v-A(cdi)v-A respectively 2 . Our first task 
is to write the operators corresponding to Eq. (J2U)) in the heavy quark effective 
theory 

For the charged modes this was already done in Ref. The operators Q\ 2 can 
be written in terms of their heavy (A, B) and light (a, b) degrees of freedom in the 
more general form 

QabAB = _ 7s)aj B 7 M (1 _ ^ (21) 



As pointed out by Grinstein and Lebed |5(, this form illustrates that there are 
symmetry relations that would allow one to extract the coupling of this effective 
Lagrangian from the measurements of B — B mixing. 

For the charged transitions this operator has to destroy a heavy and light quarks 
of flavor (A, a) = (b, u) and create a heavy and light quarks of flavor (B, b) = (c, di). 
This operator transforms as a (6^, 1#) under the chiral symmetry and, as shown in 
Ref. , Eq. (|2T|) matches in the symmetry limit of the effective theory onto 



O abAB = (3 w Tr [(6^)7^(1 - 75 )] Tr [^ bd H^(l - 7s )] (22) 

As pointed out in Ref. this same operator with (A, a) = (b, d) and (B, b) = 
(b, d) is responsible for B — B mixing so that in principle the coefficient (3y/ can be 
extracted from experiment (for this case there is an additional color factor of 8/3). 
This heavy quark symmetry relation is valid in the GL limit where — << 
Aqcd and the four velocity of the B and D* mesons is the same. For our numerical 
estimates we will use 

Pw = ^fsfoVm^B ~ (0.034 ± 0.009) GeV 3 , (23) 



where we used the decay constants = (191 ± 27) MeV |32j and f& = (225^^ ± 
21) MeV [i^l . This last one is in good agreement with the recent CLEO measure- 
ment f D = (222.6 ± 16.71^ MeV [Q. 



For our numerical estimates we will use the tree-level coefficients C\ = 0, Ci = 1. QCD correc- 
tions are known, but this difference is much smaller than other uncertainties in our calculation. 
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For the radiative decay modes that we consider in this paper there are no light 



pseudoscalars involved so we set £ = 1 in Eq. (jzzj) to obtain the matrix elements: 

<D:(v)\o abAB \B:( V ')> = wm-tt, 

<D v \O abAB \B v > = 4[3 W , 
< D*(r))\O abAB \B v > = -APwTi-v, 



<b*m\o 



abAB 



D v > 



-A(3 W 7] ■ v, 



(24) 



where rj, rj' are the vector meson polarization vectors and v the B meson velocity. 
The last two terms vanish when the condition v ■ rj = is used, we retain them 
because we will use them beyond leading order in the heavy quark expansion later 
on, as they contribute known terms of order (m& — m c )/m c . 

The neutral modes involving the weak transition B® — > D° are mediated by the 
operators Q" 2 in Eq. (J2(Jj) . To construct a matching operator in the effective theory 
we notice that they transform as (8^, 1r) under the chiral symmetry. We also need 
to extract the part of the operators responsible for destroying a heavy anti-quark 
(of flavor b) and creating a heavy anti-quark (of flavor c). A possible match for 
of the current-current form is 



= Tr 



H^l-riH® (£ f ^£) 



(25) 



However, this operator does not contribute to the processes without light pseudo- 
scalars that we are discussing. 

A possible operator that does contributes to B® — > D° transitions at tree level 
is of the form 



O 



PwTr (#fdj7,(l-75) Tr (^^(1 - 75) 



(26) 



This operator leads to matrix elements analogous to those in Eq.(j2H) with f3 w — > f3' w . 
As was the case with the charged modes, there are other possible matches but they 
lead to the same result j^. We have not found a way to determine f3' w from 
symmetry relations and must resort to the factorization model, 



P'w = -^fBfD^/m B m D = 



(27) 



An additional factor of 1/N C = 1/3 relative to (3\y occurs because the contribution 
of Q2 to these weak transitions in factorization requires a Fierz transformation and 
color rearrangement. 

For weak transitions involving the positive parity states, B — > D**, we have 
two new operators for the S doublet. The operators Q\ and Q% in Eq. ([2~U]) for the 
charged and neutral modes can be matched in the factorization approximation into 
the operators 



Q2 
Q2 



-PwTr 



(#f4)7,(l-75) Tr y(l-7fe)(fcS}°') 



PwTr (#f£jj7„(l - 75) Tr y(l - 75X6^) 



:28) 
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for transitions of the form (0 _ , — > (0 + , l + ) c . The extra minus sign is chosen 
so that the coefficients j3w and Pw are the same as those in Eqs. (J23|) and (J2"7)) if 
we take the decay constants and masses of the two doublets to be the same. This 
is approximately true for the decay constants in the analysis of Ref. 2^] based on 
QCD sum rules. 

The T = (1 + , 2 + ) multiplet does not participate in the weak transitions in this 
approximation since its decay constant vanishes |35| . 

The matrix elements obtained from Eq. (}2*Hj) (and corresponding ones for 
(0~, l~)b (0 + , l + ) c transitions) that do not involve light mesons are then 

< D*M\QI\B 1v (ti') > = -4Pww\ 

< D v \Ql\B 0v > = —4:/3 w , 

< D* v (r])\Q c 2 \B 0v > = A/3 W T]-v, 

< D v \Q c 2 \B lv (r)') > = Afarf-v, (29) 

and the same expressions with (3w — > P'w f° r the neutral modes. Finally we will 
need weak transitions from (0 + , 1 + );, — > (0 + , l + ) c . They follow from an effective 
Lagrangian like the one in Eq. (|22|) with S fields replacing H fields and produce 
matrix elements with the same sign as those in Eq. 



IV. B -» £>*7 

We now turn our attention to the single radiative decay for exclusive channels. 
Both the charged B + — > D* +/ -f and neutral B° — > D* 0/ y modes have been estimated 
before. We begin with the charged mode discussion of Ref. 0]. With the ingredients 
introduced in the previous section it is straight-forward to compute the amplitude 
from the two diagrams in Figure El For the process B+ — > D*f(r))~f(q, e) we find in 




FIG. 3: Pole diagrams responsible for B — > D*j at leading order in heavy quark and chiral 
theories. 

terms of the notation defined in the appendix, a leading order amplitude containing 
only a magnetic form factor: Fm = (Bm) LO r ^ with, 



(F M ) LO {B + - ^ + 7) = V2G F V: b V cdi f3 w e ( 



m B + — m D + 



r+=(l + ^ 2mB+ ) . (30) 
fi D + m B + + m D * 
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Similarly, for the neutral modes Fm = (-^A/)Lo r u with 



(F M ) LO (B° ^ D^) = V2G F V; b V udt P' w e( ^ ) 



r° = ( 1 + ^ , ) • (31) 



Hbo 2m bo 



I^d° rriBO + m D *o 



Using the leading order magnetic moments and the GL limit to evaluate these 
amplitudes yields r+ = —1, r° = 1/2. Combining this with exact kinematics for the 
phase space one obtains 3 

r{B+^D?i) = ^\v: b V cd fa(p w Pr+) 2 ^(mB-m D *)(m B + m D *f 

T(B^D^) = ^f\V; b V ud f a M^K-m^lK + m..) 3 . (32) 

We wish to emphasize that these results exhibit a high sensitivity to non-leading 
contributions to the magnetic moments. With exact kinematics (i.e. takin g tt ir ^ 
mr>) and using the three models for magnetic moments in Table 8 of Ref. |l7[ one 
finds 

X ~ loop 
VMD 

RQM : r+ = -4.7, r° = 0.50. (33) 

These numbers lead to predicted rates that are larger than the leading order pre- 
diction by factors between 16 and 36 for the charged modes. For the neutral modes 
the effect is more modest. Allowing the magnetic moments to vary in the ranges 
predicted by these models results in the rates 





= -6.3, 


r° 
/< 


= 0.28, 


1 /' 


= -4.7, 


r° 
1' 


= 0.38, 


' /' 


= -4.7, 


r° 
/< 


= 0.50. 



i-7 

-8 

-9 



0.7 x 10~ 7 < B(B° -> D*°-f) < 4.6 x 10' 

1.5 x 10~ 8 < B(B° S -> D*° 7 ) < 3.0 x 10" 

3.6 x 10~ 9 < B(B + -> D*+ 7 ) < 4.9 x 10" 
0.7 x 10~ 7 < B(B + -> L>*+ 7 ) < 1.0 x 10~ 7 . (34) 

There is also a significant uncertainty from the use of <?m(0) for the electromagnetic 
transitions, when in these reactions one should use gu{k 2 ~ — F"i ~ —(2 GeV) 2 ). 

We now consider two types of higher order corrections to these results that are 
counterpart to additional terms present in the pole model calculation of the neutral 
modes in Ref. jfJEij]. They correspond to the diagrams shown in Figure |H The first 



3 Our rate for the charged process is a factor of 9 smaller than the result obtained in Ref. |5| which 
corresponds to this limit. It appears that Ref. jfj incorrectly used the light antiquark charge to 
calculate the magnetic moments instead of the light quark charge. 
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D* 
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B 



D, 



b) 



B D 



D* 



c) 



B ) B i D* 
H 



d) 



FIG. 4: Non-leading contributions to B — > D*j: &)P* — P weak transitions that vanish 
in the mg = m£> limit; b-d) additional pole contributions from positive parity states. 

diagram (Figure EK) contributes beyond leading order in the heavy quark symmetry 
when we allow for irig ^ mo- It generates an electric amplitude from corrections to 
the leading order 1~ propagator fixed by reparametrization invariance 

F E (B + -> D* + j) = -V2G F V: b V cdi (3 w e ( ^ ) ( 1 - ™jf) . (35) 

\m B + - m D *+ ) V m D/ 

The next two diagrams (Figure Eb,c) involve an intermediate positive parity state 
that can be either the Dq or the D\ from the 5" multiplet. They generate an electric 
and magnetic form factor respectively that can be written as 

F E (B+ - D^l) 

Finally the last diagram (Figure HJl) involving an intermediate B\ meson contributes 
to the electric form factor as 

(Off! \ 
2 B 2 • (37) 
m B ~ m D* J 

Analogous results are obtained for the neutral mode with the obvious replacement 
V* b Vcdiflw VcbVu^Pw an d magnetic moments and masses appropriate for neutral 



= -V2G F V: b V cdi (3 w e ( ^ ) 

\m B + - m D J 

= -V2G F V: b V cd ^ w e (—&—) (l-^-). (36) 

\m B + - m D J V mjj J 
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mesons. Combining all these partial results we finally obtain 

F M {B + - A* + 7) = V2G F V: b V cd% p w e ( ^ ] (3? 

\m,B+ - m D *+J 



2 2 

2771^+ m B ~ m p+ m B - m D *+ 

fi D +m B ++m D *+ ji D + m 2 D+ m B - m D + 
,S 



F E (B + L>* + 7) = -^C^^e ( ^- 

\m B + — m B *-t 

rriB — m D *+ Hd*+ m B ~ m |>*+ A i s+ 2m B 



m B -m D + n s D+ m 2 Dt+ ji s D+ m B + m D , 



With the range of magnetic moments provided by the three models, we find 

1.2 x 1(T 5 < B(B° -> L>*° 7 ) < 3.1 x 10~ 5 , 

0.7 x 1(T 6 < B(B° S -> L>*° 7 ) < 1.7 x 10~ 6 , 

0.6 x 10~ 7 < B(B + -> D* + 7 ) < 1.0 x 10~ 7 , 

0.6 x 10~ 6 < 5(5+ -> L>* + 7 ) < 1.4 x 10~ 6 . (39) 

These ranges indicate only the uncertainty in the magnetic moments, in particular 
they do not include the uncertainty in (3w or any other parameters. For comparison, 
Grinstein and Lebed obtained B(B + — > -D* + 7 ) = 2 x 10~ 8 Q (when we correct their 
number for the missing factor 1/9). As mentioned before, the charge mode exhibits 
a large sensitivity to the value of the magnetic moments due to a partial cancellation 
between the two terms in Eq. (|30|). This sensitivity is milder when the additional 
terms are included as can be seen from the range in Eq. (}39|) . Similarly we can 
compare our result to that of Ref. B(B° — > D* 0r y) = 9.2 x 10~ 7 . Again this result 
has the same order of magnitude as the leading order contribution, Eq. (|3ip. and the 
larger number in Eq. (}3"9~j) is due to the contributions of the positive parity states. 
The prediction for B(B° — » D*°-f) in Eq. (13911 i s in fact close to the experimental 
upper bound B(B° — > _D*°7) < 2.5 x 10~ 5 j3J3] and part of the range is already 
excluded. For comparison, a recent calculation finds B(B° — > D*°^y) ~ 1.6 x 10~ 6 
pjil using a different framework. Although smaller than our range in Eq. (|39|). this 
result is not incompatible with ours given the large uncertainty illustrated by the 
differences between Eqs. (pH} and (JHSJ- 



V. 5 -» D77 AND HQET 



We are now in a position to estimate the B — > D77 amplitudes. For the process 
B — > .D77 there are 5 diagrams involving only the H doublet, shown schematically 
on FigJSJ Diagrams FigEl (d) and (e) vanish at leading order due to the condition 
Vfiig' 11 ' — v^v u ) = at the weak vertex and we left with three LO diagrams. Two 
of these diagrams have a B — D weak transition and the third one a B* — D* weak 
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FIG. 5: Pole diagrams responsible for B — > D*77 at leading order in heavy quark and 
chiral theories. 

transition. In all cases the two photons are emitted from magnetic dipole couplings 
on the external legs. 

A straightforward calculation then yields the desired amplitude for B° — > D°77. 
In terms of the form factors defined in the Appendix we find 



A — —C n 



Am y Am-E 1 Am - E 2 E X E 2 



1 

+ ^B^Di^TTT FTT + 



'E 1 (Am - Ei) E 2 (Am - E 2 ) 

B = ~, (40) 

where E\ and E 2 are photons energies, Am = m% — mo, and we have defined 

C n = 2v / 27rG F « em /3( y ^KdMl (41) 

Using these lowest order form factors we find a large range for the double ra- 
diative decay rates depending on the model used for the magnetic moments, for the 
mode with most favored CKM angles, 

0.3 x 10" 10 < B(B° -> L>°7 7 ) < 3.5 x 10~ 10 . (42) 

This large sensitivity to the input parameters is due in part to a cancellation between 
the terms involving the magnetic moments for the charm and bottom mesons. To 
illustrate this we write 

[J>d = -r^ B , (43) 
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and show in Figure |U]r(.B — > D°'y~f) as a function of r M . We use the RQM magnetic 
moment for the B and we normalize the rate to its value when r M ~ 2.87, the value 
of no in the RQM. 

25 n — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — . — ■ — ■ — ■ — ■ — n 



20 - 
_ 15 : 

a. 
U 

^ 10 : 

5 : 




1 2 3 4 5 

FIG. 6: T(B° — > Z)°77) as a function of r^. We use the RQM magnetic moment for the 
B and we normalize the rate to its value when ~ 2.87, the value of [id in the RQM. 



B {H,S} ( {D 1; D*} < D 
X ' 1 



a) 



B ( { B lt B*} ( { H,S } D 



b) 




FIG. 7: Additional diagrams involving particles from the (0 + , 1 + ) and (l + ,2 + ) doublets. 
H(S) stands for either of the H(S)-multiplet members 

We next consider two other types of contribution. First the non-leading contri- 
butions from diagrams d and e. They arise from corrections to the 1~ propagator 
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B S ( D] ( D 
X ^ i 



B < Bi < S D B(B|,B* Di,D*<D 
3 '- X — J X ■> 



a) 



b) 



c) 




FIG. 8: Additional diagrams involving particles from the (0 + ,l + ) and(l + ,2 + ) doublets. 
S stands for either of the S-multiplet members 



that are fixed by reparametrization invariance resulting in 

1 1 



n — ^ n 



HdUd* m% - m 2 D 



Am 



in 



D 



+ 



Am — Ei Am — E 2 



HbUb* rn 2 B - m 2 D 



E X E 2 



in 



B 



(44) 



D 



C n 

T 



HdUd* m 2 B -m 2 D ( 



Am 



in 



D 



Am — E\ Am — E 2 



jj, B ^B*(E 2 - Ex) m 2 B - m 2 D 
E 1 E 2 {E l + E 2 ) m% 



Additional contributions arise when the positive parity states appear as in Figure 
The contributions from the S multiplet are 



.4 



B 



fi s D fi D m| - m 2 D 



Am 



hUHd + Ho- 



rn 



i) 



Am — Ei Am — E 2 
1 1 



HbHb m B - m 2 D 



E\E 2 



rn 



B 



Am 



D 



Am 

C n r 

2 [ Am 

i4(Eh - go 

E\E 2 (Ei + E 2 



Hd /„ s m 2 B -m 

[flu - /i D 



Am — Ei Am — E 2 
1 1 



q / q mi,— mi-, • 

Hb(j*b- Hb^^, 
E\E 2 



in 



D 

m 



+ 



D - 



m D 

(Hb + Hb- 



Am — Ei Am — E 2 
1 1 



Am — Ei Am — E 2 



m 



D 



III 



B 



(45) 



S 2 2 

Hb i , , S 171 B ~ 171 D - 



Ei E^ 



m. 
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The contributions from intermediate T doublet states are 



B = AC n 
C = -2Q 
D = -2G 



Kr'D) 1 


1 


Am \ 


Am - E x 




m 2 B — m 2 D 


Am 


m 2 D 


>S) 2 


m 2 B — m 2 D 


Am 


m 2 D 



1 



Am - E 2 
1 

Am - E x 
1 



Ml! 

EiE 2 
1 



(46) 



Am - E 2 
1 



+ 



« 2 



m B — m 



D 



m 



B 



Am — Ei Am — E 2 



E\Ei 



(^ B ) 2 (E 2 -Ei)m 2 B 



m 



D 



E 1 E 2 (E 1 + E 2 ) 



m 



B 



Additional contributions from SS and TS multiplet transitions as in Figure |H1 
We split them into those from diagrams with two members of the S doublet: 



B 



Cn 



+ /4/4( 



Am 
1 



( 



+ 



fi B fiB* 



Am — Ei Am 
1 

+ 



E 2 



E X E 2 



q _ C n M S D [lp 
" 2 1 Am 



Ei(Am-Ei) E 2 (Am-E 2 
1 1 



Am — Ei Am — E 2 



+ 



E\E 2 



(fi s B fi D - fi s D fi B )( 



Ei(Am-Ei) E 2 (Am-E 2 



■)]• 



(47) 



D 



2 ^ Am \Am-Ei Am - E 2 
+ (jt&D + 1* S d Hb){ 



H b Hb(E 2 - Ei 
EiE 2 (Ei + E 2 ) 

-)], 



(Am - Ei) E 2 (Am - E 2 

and those from diagrams with one member of the S doublet and one member of the 
T doublet: 



D 
C 
D 



2C n 

a/3 L Am 
V3 



T TS 2 

HdUd m B 



m 



D 



Am 



m 



v 



+ 



1 



Am — E\ Am — E 2 
1 



T TS 2 
f^B^B m B 



m 



D 



EiE 2 



m 



B 



Am - Ei 
1 



Am - E 2 
1 



Am V Am — Ei Am — E 2 



+ 



Ei E 2 _ 

H T B fi T B S (E 2 - Ei) 
EiE 2 {Ei + E 2 ) 



(48) 



In a similar manner we obtain the result for — > D ±r y^ with the replacements 
V cb V ud f3' w -> V ub V cd [3 w , for B° s D° 77 with V ud V us and B ± £±77 with 

VcbVudP'w ~* VubVcsPw- 

Numerically we find the following ranges for the branching ratios when we vary 
the magnetic moments over the range predicted in the three different models 

1.7 x 1(T 8 < B(B° -> D°7 7 ) < 8.0 x KT 8 , 

0.8 x 10" 9 < B(B° S -> L> 77 ) < 4.3 x 10~ 9 , 

0.6 ■ 10~ n < B(B ± -> £±77) < 2.0 ■ 10" 11 , 

1.1 ■ 10~ 10 < B(B ± -> Dfjj) ~ 3.6 ■ 10~ 10 . (49) 
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Once again these ranges include only variations of the rates with the magnetic 
moments in the %LM, VMD and RQM. They do not include other uncertainties 
such as that in the value of @w 

It is instructive to examine some of the features of the differential decay rates 
and we do so for the mode with most favorable CKM angles, B° — > D°77. We 
first plot in Figure 01 the normalized differential decay rate as a function of z, the 
dimensionless photon pair invariant mass defined in Eq. (|A7[) . The distribution does 
not have any thresholds as our calculation does not include absorptive parts. It is 
peaked at the higher invariant masses. To evaluate the limits of the heavy quark 
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FIG. 9: Normalized differential decay rate for B — > D 77 as a function of the photon 
pair invariant mass. 

expansion we first show in Figure ITU1 the double differential decay rate as a function 
of the energies of the two photons as a density plot. The darker regions correspond 
to the most populated ones. The distribution is dominated by the region in which 
both photons tend to have similar energy between 1 and 2 GeV. This indicates that 
the approximation of constant magnetic moments for the photon emission vertices 
is slightly better for the double radiative decay modes than it was for the single 
radiative decay mode where E y ~ 2.3 GeV. However, it is clear that a substantial 
uncertainty remains due to this approximation. 

Finally we show again the normalized differential decay rate as a function of 
uj = v ■ v' in Figure Here we see that the distribution is peaked at uj ~ 1.1, 
not too far from the symmetry limit. The heavy quark expansion should be better 
behaved for these modes than it is for the single radiative decays where o> ~ 1.5. 
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FIG. 10: Density plot for the double differential decay rate for B° —>■ D°j"f as a function 
of the two photon energies. 

A. Additional Contributions 

We turn our attention to potential contributions with a different topology that 
have not appeared so far. Specifically we have in mind contributions in which the 
weak decay B — > DX is followed by the electromagnetic X — > 77 vertex. We do not 
have a systematic way to include these contributions, but we illustrate them with 
a few examples. In any case, if there is a dominant contribution of this form, the 
photon pair invariant mass would be concentrated around M x and would be easy 
to isolate experimentally. 

An example of this topology with a cc resonance is B — > Drj c — > .D77. Using the 
narrow width approximation one finds T(B — > -D77) = Y[B — > rj c D) ■ B(r] c — > 77). 
Experimentally it is known that B{rj c — > 77) = 4.3 ■ 10 -4 0. For the weak vertex 
we can use a recent pQCD calculation as an illustration, B(B — > rj c D) = 1.28 • 
10~ 5 This then leads to a contribution B(B° —> D°77) ~ 5.5 • 10~ 9 smaller 

than our result in Eq. H2 A different estimate B(B r] c D) = 1.52 • 10" 7 
would make this contribution even smaller. Another charmonium resonance with a 
measured two photon width is the Xco(l-P)(0 + ). Since the measured branching ratio 
is i3(xco(l-P)(0 + ) — > 77) = 2.6 • 10~ 4 , similar to that for r\ c — > 77 we expect at most 
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FIG. 11: Normalized differential decay rate for B° — > _D°77 as a function of a> = i> • i/. 

a similar contribution to the double radiative S decay. 

We expect much larger contributions from the light pseudoscalars. This can be 
illustrated using the measured rates for B(B° — > D°X), X = tt , 77, 77' which are all 
at the 10 -4 level, and the respective two photon widths which range from nearly 
100% for the ir° to a few percent for the rj'. These contributions would however, be 
easily separated experimentally as the two photon invariant mass distribution would 
correspond to sharp peaks at the respective mx- 

VI. SUMMARY AND CONCLUSIONS 

We studied the weak radiative decays of B mesons that occur at tree-level in 
the Standard Model. We presented a numerical estimate for the inclusive double 
radiative decay b — > X c 77 based on a free quark decay calculation. Our estimate 
indicates that this mode is about an order of magnitude larger than the double 
radiative penguin mode b — > ^77. As such it can be studied at future Super-B 
factories in experiments designed to observe the double radiative penguin mode, and 
in fact will be a significant background to that mode. 

We reviewed the HQET formalism in connection with the study of the single 
radiative exclusive modes of the form B — > D*^j. We first reproduced existing results 
in the literature for both charged and neutral modes with a leading order calculation 
including only the (0~, 1~) doublet as intermediate states. These calculations exhibit 
a large sensitivity to the value of the electromagnetic couplings due to a partial 
cancellation. We then improved these lowest order results by including certain 
known terms of order (m^ — m c )/m c as well as by introducing the positive parity 
doublets S and T as intermediate states. These two ingredients significantly enhance 
the predictions as they remove the cancellation that occurs at lowest order. The 
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framework is a good approximation in the m& ~ tti c limit, but significant corrections 
are expected for the physical values of quark masses. 

Finally we extended the calculations to the double radiative decay modes of the 
form B — > .D77. Once again we found significant cancellations between the lowest 
order terms and a much larger rate when the positive parity states are included in 
the calculation. The mode with the most favorable CKM angles is predicted at the 
level, comparable to predictions for B — > .K77. 

We expect the HQET formalism to work best in the case when the velocity of the 
heavy hadrons remains constant during the transition. Whereas this is kinematically 
impossible in the single radiative decay modes, there are regions of phase space in 
the double radiative decay modes where this could be tested (in principle at least). 
We illustrate these regions with plots of the relevant differential decay rates. 
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APPENDIX A: KINEMATICS 



The most general amplitude for decays of the form M — > V(r))^y(q, e) consistent 
with electromagnetic gauge invariance can be written in terms of two form factors. 
Labeling the momentum of the initial state M with its four velocity in its rest frame, 
p = Mmv, and denoting the polarization of V, 7 by rj, e respectively 

M(M^ Vj) = e*" [iF M e^v v q a ri^ + F E (q ■ rfv^ - q ■ vrQ] . (Al) 

Summing over the photon and vector meson polarizations the partial decay rate is 
given by 

m - ^7) = (\ F m? + \F E \ 2 ) , (A2) 

For the heavy meson formalism in which the meson fields are normalized as 

< M(v',k')\M(v,k) >=2v°5 vv ,(27i) 3 5 3 (k-k'). (A3) 
the decay rate becomes instead 

r(M^ 7 ) = gi|( W + W ), (A4) 

where we have used exact kinematics, otherwise in the heavy quark symmetry limit 
Ey = Mm as well. 
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The amplitude for decays of the type M — > M /r y~f, with M, M' pseudoscalar 
mesons, 



M{M(p) - M'ipMh, e 1 ) 1 {k 2 , e 2 )) = e^(A; 1 )e* l/ (A; 2 )M^(p, h, k 2 ) 



(AS) 



is well known from the kaon literature [41|. The most general decay amplitude 
M tM1/ (p, ki, k 2 ) consistent with electromagnetic gauge invariance and Bose symmetry 
contains four form factors. We write them here in terms of the velocity of M in its 
rest frame, so that p = M^v as it appears within the HQET. 



ill 



/<!/ 



M 2 M 



{k 2fl k lu - g^kx ■ k 2 ) 



. C(z,y ) 
1 Ml 



k\ k 2 



(A6) 



+ % 



—j^i ^ \ v ' hvvfa^ + v ■ k 2 v IJl k lv - v ■ kiv ■ hg^ - k x ■ k 2 v^v u ] 
D(z,y) 



M 2 M 



v ■ he^apk^v 13 + v ■ k 2 e IJVa ph^v p + (v M e^ 7 + v^e^^k^k^v 1 



where 



v • (Aii - k 2 ) (Aii + k 2 )' 

y = — , z - 



M M 1 ~ M' M 1 M M (A?) 

The relation between these dimensionless variables and the energy of the two photons 



is 



Ex = ^M M ((z + 2y + l)-r 2 ) 
E 2 = l -M M ((z-2y+l)-r 2 ) 



(A8) 



Recently, Hiller and Safir |42j have claimed in the context of B — > K'yy that there 
are three additional form factors, 



.C+(z,y) 

Ml 
C-(z,y) 



+ i 



Ml 



+ i 



D-(z,y) 



Ml 



k\ ■ k 2 e^ vaf3 {kx + k^)v p + (k 2/M e va p 1 + k^e^^k^v 1 (A9) 
k\ ■ k 2 e^ ua p{ki k 2 )v^ {k 2fl € ua /3^ ki u e lia p^)kik 2 v' Y 



It is well known, however, that these form factors are not independent and can be 
reduced to the ones in Eq. (|A6[) . This follows from the fact that in four dimensions 
there are at most four linearly independent four vectors, and this gives rise to the 
Schouten identity [I^f. In this case all three terms in Eq. (|A9|) reduce to the second 
form factor in Eq. (|A6J) so that 4 

= (m mV C + {z, y) + ^i(l + z -r 2 ) C~{z, y) - D~(z, y)j c^AfA&AlO) 



The authors of Ref.|42| agree with us and have added a note to this effect in their paper. 
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The physical region in the dimensionless variables z and y is given by 

0<\y\< ^(l, r \z), < z < (1 - r) 2 , (All) 

with 

A(a, 6, c) = a 2 + b 2 + c 2 - 2(ab + ac + be). (A12) 

Note that the invariant amplitudes A(z, y), B(z,y) and C(z,y) have to be sym- 
metric under the interchange of k\ and &2 as required by Bose symmetry, while 
D(z,y) is antisymmetric. Using the definitions (|A6|) the double differential rate for 
unpolarized photons and conventionally normalized meson fields is given by (in the 
rest frame of M) 

d " T Mm [z 2 ( \A + B\ 2 + \C\ 2 ) (A13) 



dydz 2 9 tt 3 

+ (\B\ 2 + \D\ 2 ) (y 2 -\\(l,r 2 ,z)n 

With HQET normalization for the meson fields and with the two photons re- 
taining their usual normalization Eq. (|A13|) is replaced by 

d 2 T _ M 2 M E M , [z2{ lA + Bl 2 + lcn (AM) 



dydz 2 9 tt 3 

+ (|5| 2 + |^| 2 ) (y 2 -\x(l,r 2 ,z)n 

of course, in the heavy quark limit Em> — * M' as well. 
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